Tunneling- induced renormalization in interacting quantum dots 
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We analyze tunneling-induced quantum fluctuations in a single-level quantum dot with arbitrar- 
ily strong onsite Coulomb interaction, generating cotunneling processes and renormalizing system 
parameters. For a perturbative analysis of these quantum fluctuations, we remove off-shell parts of 
the Hamiltonian via a canonical transformation. We find that the tunnel couplings for the transi- 
tions connecting empty and single occupation respectively single and double occupation of the dot 
renormalize with the same magnitude but with opposite signs. This has an important impact on 
the shape of the renormalization extracted for example from the conductance. Finally, we verify 
the compatibility of our results with a systematic second-order perturbation expansion of the linear 
conductance performed within a diagrammatic real-time approach. 

PACS numbers: 73.23.-b,73.23.Hk 



I. INTRODUCTION 

Few-electron quantum dots are paradigmatic systems 
to investigate the effects of Coulomb interaction and 
quantum fluctuation in nanoscalc systems^— The theo- 
rist's workhorse in this field is the so-called Anderson- 
impurity model, consisting of a single spin-degenerate 
level with onsite Coulomb repulsion U, tunnel coupled 
to non-interacting leads. In this type of systems, the 
energy scale associated with Coulomb repulsion is usu- 
ally large and consequently interaction effects cannot be 
treated within a perturbative scheme. However, when 
the tunnel coupling between dot and leads is smaller 
than temperature or voltage, a perturbation theory in 
the tunnel coupling strength, here denoted by T, can be 
successfully employed. Transport in lowest order (some- 
times referred to as the sequential tunneling) can in 
many cases be understood by a straightforward master- 
equation approach with transition rates computed by 
means of Fermi's golden rule^ see for example Ref. @- 
In this so-called orthodox theory, energy conservation al- 
lows only for transitions between energetically degener- 
ate states, and we refer to them as classical or on-shell 
transitions. 

Higher-order transport corrections are associated with 
quantum fluctuations or off-shell transitions. In second 
order in T, there are two different types of quantum- 
fluctuation corrections. First, there are cotunneling pro- 
cesses. They consist of a sequence of two tunneling events 
with an intermediate virtual state (i.e., the energy of 
the intermediate state is not equal to the energy of the 
initial and final state). Cotunneling dominates trans- 
port in the Coulomb-blockaded regions^ Some of the 
cotunneling processes, e.g. spin flips, occur already at 
zero voltage, other effects were predicted to contribute 
at large bias voltage, e.g. inelastic cotunneling^ and pair 
tunneling £t~— The second type of quantum- fluctuation 



corrections are described by renormalizations of the sys- 
tem parameters, namely the dot energies and the tunnel- 
coupling strengths. They lead to shifted positions and 
modified tunneling strengths at the transport resonances 
as a function of gate voltage. In many cases, these renor- 
malizations yield only small corrections to the lowest- 
order contribution to, e.g., the conductance. There are, 
however, also scenarios in which the renormalization cor- 
rections are of crucial importance, because the lowest- 
order contribution either vanishes (as it is, e.g., the case 
for certain schemes in adiabatic pumping through single- 
level quantum dot o 13 ' 14 ) or provides a flat background 
only (as, e.g., in the maximal linear conductance through 
a metallic single-electron transistor, where a logarithmi- 
cally temperature-dependent conductance indicated mul- 
tichannel Kondo behavio r 15 ' 16 ). 

Any theoretical approach that takes systematically all 
contributions to a given order in perturbation expansion 
into account will automatically include both cotunnel- 
ing and renormalization effects, and, thus, go beyond 
pure cotunneling as obtained from standard second-order 
perturbation thcoryj&i£ This is true, e .g.. for the dia- 
grammatic real-time approach, see Refs.ll8l that was de- 
veloped to describe transport through quantum-dot sys- 
tems. In this case, however, the renormalizations can 
only be read off a posteriori by bringing the final re- 
sult of the considered quantity (e.g., the linear conduc- 
tance) into a form that allows for an identification of the 
renormalizations. We show in this paper that the lin- 
ear conductance, as one possible measurable quantity to 
be considered, reveals these renormalization effects; how- 
ever, care has to be taken, since an identification relying 
only on the conductance may not be unique. Therefore, 
we introduce in this paper an a priori procedure to clas- 
sify and evaluate the various quantum- fluctuation correc- 
tions. This procedure is based on a canonical transfor- 
mation that removes off-shell parts of the Hamiltonian 
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and, simultaneously, generates new transitions describ- 
ing cotunneling as well as renormalization of the system 
parameters. The main emphasis of this paper is on the 
derivation of the explicit expression for the renormaliza- 
tion of the tunnel coupling strength. 



II. MODEL AND METHOD 

As a specific example, we consider a single-level quan- 
tum dot with Coulomb interaction weakly coupled to two 
reservoirs. Since we are interested in the linear conduc- 
tance, a small bias voltage can be applied between the 
two reservoirs. The system consisting of quantum dot 
and leads is described by the Hamiltonian 

H = if dot + -^tunnel + H\ ca( i . (1) 

The single-particle level spacing in the dot is assumed 
to be larger than any other energy scale (temperature, 
Coulomb interaction, transport voltage) such that only 
one energy level needs to be taken into account. Hence, 
the dot can be described by the single-level Anderson 
model 

H dot = ed Ua + Un^n x . (2) 

The creation (annihilation) operator for an electron 
with spin a on the dot is given by $ a {d a ), and n a = 
d\,d a is the corresponding number operator. The on- 
site repulsion U (as found from the constant interac- 
tion model^) describes the energy cost for double oc- 
cupation and stems from Coulomb interaction. Tun- 
neling of electrons between dot and leads is taken into 

account by Jftunnei = Y, a ,k,a V a cj L,k,a d a + n - c - We as- 
sume a momentum- and spin-independent tunnel ma- 
trix element V a and define the creation (annihilation) 
operators <y k cr (c a krr ) for electrons with spin a and mo- 
mentum k in lead a = L, R. The leads' Hamiltonian 
is given by ff lcads = E a fe, CT e a,fe4 ifei0 .c a>fe . The chemi- 
cal potential of the two leads differs by the applied bias 
— (J>r = —eV, with e > being the electron charge. 
We assume that the density of states p a in the leads is 
constant in the window relevant for transport and define 
the tunnel coupling strength r a as r Q = 2irp a \V a \ 2 and 

r = r L + r R . 

The (reduced) Hilbcrt space for the quantum dot is 
spanned by the states |0) for an empty dot, \a) for a 
singly-occupied dot with spin a =t, 4-, arL d |d) for a 
doubly-occupied dot. The corresponding energies are Eq, 
E a , and E^. The (high-dimensional) Hilbert space of the 
full problem, on the other hand, is spanned by the many- 
body eigenstatcs \n) of the dot decoupled from the leads 
and of the leads, with energy E n (containing both the 
energies of the lead and the dot electrons). 



III. CANONICAL TRANSFORMATION 

For a systematic analysis of quantum-fluctuation ef- 
fects due to tunneling, we split the Hamiltonian into three 
parts; a term in the absence of tunnel coupling, Hq, and 
two different types of tunneling, H\ and H%, 

H = H + H 1 +H 2 . (3) 

The dot and the reservoirs in the absence of tunneling 
are described by 

H = Y,E n \n){n\. (4) 

n 

The tunneling part of the Hamiltonian introduces cou- 
plings between different eigenstates \n) and \n'). We dis- 
tinguish on-shell contributions, 

Hx = Y,Vn'n\n')(n\5 EnjEn , , (5) 

nn' 

which couple states of the same energy, E n = E n i, from 
the off-shell parts, 

H2 = J2 V n>n\n'){n\(l-5 En , En ,) , (6) 

nn' 

which connect states with different energies, E n ^ E n >. 
In both cases, we have used the abbreviation 

V n 'n = (n'\H tunncl \n) . (7) 

The rates for classical (on-shell) transitions between 
two states |n) and |n') are obtained via Fermi's golden 
rule in first order in the tunnel coupling T corresponding 
to second order terms in V, 

Wn'n = |K'n| 2 5 En ,E n , , (8) 

n 

where the Kronccker delta ensures energy conservation. 
As a consequence, only the transitions described by the 
on-shell part Hi need to be considered. In order to find 
the transition rate between two states of the dot subsys- 
tem, an average over all possible initial lead states has 
to be performed. The off-shell part H2 contributes to 
higher orders in the tunneling only. It is this part that 
describes quantum fluctuations. 

The aim of this section is to remove Hi and account 
for its effect by renormalizing the system parameters en- 
tering the on-shell part H + Hi and by generating new 
transitions. In the following, we derive this renormal- 
ization to lowest order in T, i.e., in second order in the 
tunneling matrix elements V. To this end we perform a 
canonical transformation^— 

H = e~ iS He iS (9) 

with the Hcrmitian operator S being chosen such that Hi 
is eliminated. This is achieved by Hi+i[Ho, S] — 0, i.e., S 
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is linear in H2] the explicit form of the matrix elements 
of S is given by (m\S\n) = i(m\H 2 \n) / (E m — E n ), for 
E m 7^ E n , and is equal to zero otherwise. With this 
condition, the transformed Hamiltonian H, expanded up 
to third order in the tunnel coupling V, reads 

H = H tt + H x - l - [S, H 3 ] - i [S, Hi] - \ [S, [S, Hi]] 
~[S,[S,H 2 ]] . (10) 

We split H = Hq + Hi + H2; again, into a diagonal 
part Hq, an off-diagonal but on-shell part Hi, and an 
off-shell contribution H2 ■ The new effective model is ob- 
tained by dropping H2, i.e., by dropping —i [S,H±] and 
all other non energy-conserving contributions which ap- 
pear in Eq. (fTU)) . that would contribute to processes in 
yet higher orders only. In the next three subsections we 
analyze the effect of the different corrections contributing 
to H, 




e/r 



FIG. 1. Level renormalization Se as a function of the bare level 
position e in units of T, for zero bias voltage. We furthermore 
choose ksT = 2F (full line) respectively feT = (dashed 
line) and U = 30I\ 



A. Energy renormalization 

The diagonal term, Hq, contains the renormalized en- 
ergies E n = E n + SE n (for the combined system of dot 
plus reservoirs), with 



SE n 



V \(m\H 2 \n)\' 
^ E — E 



(ii) 



The renormalization depends on the initial state n. Since 
we treat the leads as a reservoir in equilibrium, we 
average over the reservoir part of the initial states n 
according to the Fermi distribution function f a (u>) = 
{1 +exp[/3(w — p a )]} _1 where (3 = l/(k-B,T) is the inverse 
temperature. This leads to the following renormalization 
of the dot energies, 



SE = 
SE a = 
5E d = 



la 

2tt 



e — uj 

l-/«( w ) , /a( w ) 



uj — e e + U — ui 
uj-e-U ' 



(12) 
(13) 
(14) 



Energy renormalizations of this type have been discussed 
in more complex quantum dot systems, where a non equi- 
librium of a spin or pseudo spin occupation can occur, 
see e.g. Ref. [24|. We refer to this in more detail in Sec- 
tion HTC] The integrals of Eqs. JT2]) to {H} - and the ones 
which will be subsequently discussed - are regularized 
through Cauchy's principal value (which is equivalent to 
adding +i0 + in the denominators of the integrands and 
taking the real part after integration). This ad hoc regu- 
larization procedure has been addressed rigorously in the 
context of a T-Matrix approach in Ref. [2f|. 



The addition energies are defined as e = E a — Eq and 
e + U = Ed — E a . From the expressions in Eqs. (fT2")) - 
(|14p one can extract the renormalization contribution to 
the addition energies, namely Se = 8E a — SEq and SU = 
SEd — SE a — 8e, which are given by 



Se = — — — / dui 



^ 2tt 



uj — e 



i + U-LJ 



SU = 0. 



(15) 
(16) 



The integral can be performed analytically leading to the 
energy level renormalization 



5e = J2 T * [0a(e + EO-0«(e)] 



(17) 



^Re* 



2n 



where ^(x) is the 



with c/> a (x) = 
digamma function. 

The energy renormalization as function of the level po- 
sition in the absence of a bias voltage is shown in Fig. [1] 
At zero temperature, Se = ^ In | diverges logarith- 
mically at the degeneracy points for empty and single 
occupation, e = 0, and for single and doubly occupation, 
e = —U. At finite temperature, the divergency is cut 
off. The renormalization Se vanishes at the particle-hole 
symmetry point, e = —U/2. Away from e = —U/2, the 
sign of Se is such that the effective resonance positions of 
the system are always shifted away from the particle-hole 
symmetry point. 



B. Cotunneling 

The off-diagonal, on-shell part Hi of the transformed 
Hamiltonian describes tunneling transitions. We start 
with discussing the contributions occurring in second or- 
der in the tunnel matrix element V. They stem from 
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— (i/2) [S, H 2 ] and are associated with cotunncling pro- 
cesses with tunnel amplitudes (with n n') 



(n'\H 2 \m)(m\H 2 \n) 



En 



(18) 



We remark that n and n' have the same parity in the dot- 
electron number, in contrast to the bare tunnel coupling 
V n ' n that is only non-zero if n and n' have opposite parity. 

The transition rates for these cotunneling processes 
are obtained by plugging the amplitudes into the Fermi 
golden rule, w n > n = x \5V n > n \ 6(E n < - E n ), and, after- 
wards, averaging - as we did before - over the initial 
reservoir occupation. The rates describing transitions in 
which the dot charge changes by 2e are given, in their 
integral form, by 



w d0 



2nh 



dw/+M/+,(2e + tf-a;) 



1 



1 



^ w — e u> — e — U 

v r / 



■fxl 



2nh 



duf-{u)f-,{2e + U-w) 
1 ^ 2 



oj — e ui — e — U 



(19) 



(20) 



where we used the notation /^"(w) = / Q (w) and /~ (cj) = 
1 — / q (cli) and a, a' denote the leads involved in the tun- 
neling process. Furthermore, there are rates for processes 
in which the spin of the dot is flipped via an empty or 
doubly occupied dot as intermediate state, 



r r 1 

A a 1 - a' 

2nh 



1 



1 



uj — e ui — e — U 

(21) 

Here, a denotes the spin opposite to a. There are also 
rates, that leave the state of the dot unchanged. This 
happens when an electron enters from one lead and leaves 
to possibly another one. The rates are given by 



1 00 



du> 



/ (w - e) 2 



2nh 



1 



+ 



1 



(w-e) 2 (oj-e-Uy 



„,«» 

^dd 



r r / 

2irh 



du 



(w - e - C/) 2 



(22) 

(23) 
(24) 



Again, all the integrals are regularized by adding +i0 + in 
all resolvents and taking the real part after integration. 
This integral form is convenient to trace back how many 
electrons have left or entered a specific lead a: the factor 
fa( u ) indicates that an electron has left and /~ (w) that 
an electron has entered lead a. The analytic evaluation 
of the integrals is straightforward; the result is given in 
the appendix. 



C. Tunnel-coupling renormalization 

Finally, we deal with the contributions to the off- 
diagonal, on-shell part Hi, which occur in third order 
in V. These contributions rcnormalize the tunnel matrix 
elements V n 'n- We get 



SVn'n = - 



(n'\Hi\n) 



E 



\{m\H 2 \n)\ 2 + \(ri\H a \m)\ 2 



E 



^ {E m - E n f 

n'\H 2 \m') (to' I H 2 \ m) (m\H 2 \ n) 
(E m > - E n )(E m — E n ) 



(25) 



Again, we subsequently average over the occupation of 
the leads in order to get the renormalization of the tun- 
nel matrix elements of the bare Hamiltonian, V a . Both 
the situations where an initially occupied or an initially 
empty state in one of the leads is required appear; this 
leads to contributions with the weight f^,(u>) and f~,(uj). 
After combining them and making use of the fact that 

f du!-? r -n+vi = 0, wc find for the renormalization of 

the transition amplitude connecting an empty with a 
singly-occupied dot (with dot excitation energy e) 



1 \ ^ r a > 

2 ^ 27 



duf a ' (w) 



1 



1 



oj-e-U 



(26) 

The subscript e indicates that the transition between an 
empty and a singly occupied dot with an excitation en- 
ergy e is considered here. The conservation of hermitic- 
ity of the canonical transformation yields SV* e /V* = 
5V a . t /V a . Most notably this renormalization is different 
from the transitions connecting a singly-occupied with a 
doubly-occupied dot (with dot excitation energy e + U). 
To be more specific, we find that it is equal in magnitude 
but has an opposite sign compared to the previous one, 



sv a 



e+U 



8V r 



a,e+U 



V* 



6V a ,e 



(27) 



As a consequence, the renormalization of r Q is given by 



r„ 



Y,r a ,Q a ,(e,U) , (28) 



shows that the ratio 



with Q a ,(e,U) = [<f/ a ,(e + U) 

0o.'(e)}]. The result of E q- 
(5r a £ /r Q is the same for coupling to both leads and 

therefore in the figures we will suppress the index a. 
Furthermore, in the zero-bias case this ratio does not 
depend on the asymmetry of the bare couplings but 
only on r = Tl + Tr,. The renormalization of the 
tunnel couplings <5r Qj£ and 5T aje+ u in the absence of 
a bias voltage is shown in Fig. [2] At zero tempera- 
ture, the renormalization of tunnel couplings 5T a ^ e = 



' 2tt 



e+U + 

and e 



A In 



e+U 



-U. 



displays 1/x divergencies at 
At finite temperature, these diver- 



gencies are cut off. The sign of the renormalization is 
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IV. LINEAR CONDUCTANCE 

In the following, we show that the analysis presented 
above can fully explain the quantum-fluctuation correc- 
tions to the dimensionlcss linear conductance g = Gh/e 2 
with G = (dI/dV)\v=o through a single-level quantum 
dot. By using current conservation 7l = — Jr, the cur- 
rent can be written in a symmetric form, 1 = 1^ = 
(Tr/l — Tl^r)/ (Tl + Tr), where I a is the electric current 
flowing from lead a into the dot. 



First order 



FIG. 2. (color online) Renormalization of the tunnel cou- 
pling strengths for the two transitions between empty and 
singly-occupied (ST t , black) and between singly- and doubly- 
occupied dot (5r e +[/, red) as functions of the level position 
e in units of F, for zero bias voltage. We furthermore chose 
feT = 2r (full lines) respectively IcbT = (dashed lines) and 

u = 3or. 



such that the excitation (either e or e + U) that is closer 
to the Fermi energy of the leads acquires a stronger tun- 
nel coupling if the dot is preferably singly occupied and 
a weaker one for preferred empty or double occupation. 
The sign for the renormalization of the tunnel coupling 
associated to a certain excitation (cither e or e + U) is 
opposite if the other excitation (namely e + U or e) is 
close to the Fermi energy. 



D. Comparison with poor man's scaling 

Deriving an effective low-energy Hamiltonian with 
rcnormalizcd system parameters is also the central idea 
of renormalization group approaches. Its simplest ver- 
sion, poor man's scaling, iteratively integrates out the 
high-energy degrees of freedom.] 26 ! 27 It is thus possible to 
sum up the leading logarithmic contributions. However, 
this approach does not properly account for the sublead- 
ing terms. Applied to the single-level Anderson impurity 
model for U » lei, one obtains 



e ^ 2tt n max{|e|,fc B T} 



(29) 



We start by considering the current I a expanded to 
first order in the tunnel coupling. It is obtained by eval- 
uating the first-order tunneling rates w X ' X that describe 
the change of the dot state from \ to x\ see Eq. (©, 
and multiplying them with a factor +1 ( — 1) when an 
electron enters (leaves) the dot during the transition. Fi- 
nally, these rates multiply the (zeroth-order in the tun- 
nel coupling r) probability p x {V) to find the dot in state 
X € {0,t, 4-, d}, when a bias voltage V is applied. After- 
wards, we linearize in V and arrive at the dimensionlcss 
linear conductance 
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(i) 



<T=t4 



r T ,r 



L-L R 



r L + r R 



/'(e) (Pa+Po) 



r Tj r 



Li R 



f'(e + U)( P<J +p d ) 



•(30) 



Here, p x denotes the (zeroth-order) equilibrium (V = 0) 
probability to find the dot in state %■ It is given by the 
Boltzmann factors p x = exp(— (3E X )/ ^ , exp(— flE' x ). 
The dummy index e and e + U attached to the ratio 
of the tunnel couplings indicates the transition (e for 
empty/single occupation and e + U for single/double oc- 
cupation). The tunnel couplings are the same for both 
transitions. However, as discussed above, they renormal- 
ize differently, which is the motivation for distinguishing 
them. 



B. Second order 

We subsequently discuss corrections to the linear con- 
ductance in second order in the tunnel coupling T. Ac- 
cording to the above discussion about the different con- 
tributions to the quantum- fluctuation corrections, we ex- 
press the second-order linear conductance as the sum 



which qualitatively describes the correct low-temperature 
energy renormalization. The exact determination of the 
high- and low-energy cutoff appearing in the argument of 
the logarithm is, of course, not possible within the poor 
man's scaling approach. Furthermore, it cannot address 
the tunnel-coupling renormalization at alli^l 



o (2) - o (2) 

y — fcot 



9 (2) 



E 



(2) 
5ren,r o 



(2) 
3rcn.r o 



, (31) 



containing cotunneling and different types of renormal- 
ization terms. The cotunneling part is obtained in the 
same way as first-order transport with the difference that 
cotunneling rates, see Eqs. (JTUJ) to (|2~i|) . instead of the 
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sequential-tunneling rates are used. In order to calcu- 
late the current in lead a, these rates have again to be 
weighted with factors 0, ±1, respectively ±2, depending 
on the number of particles transferred from or to lead 
a for a certain choice of a and a' . The renormalization 
parts are given by 



(2) 
9rcn,X 



dX 



5X 



(32) 



with X 6 {e, r Q!C , r aie+ [/}: where we treat T a ^ and 
T aie+ u as independent parameters. The values for 5e, 
<5r Qj£ and 5r aye+ u are taken from Eqs. (fTT)) and (|28l) . 

To verify the validity of this result, we recalculate 
the conductance within a diagrammatic real-time ap- 
proach^ The virtue of this approach is that it allows 
for a systematic perturbation expansion that does not 
require the identification and the separation of different 
quantum- fluctuation contributions. The analytical ex- 
pression that we obtain for the first- and second-order 
conductance within this alternative approach is identical 
to the one obtained from Eqs. (|3T))) and (f3"Tj) . 

The results for the conductance are plotted in Fig.[3ta) 
as a function of the level position. The black, dashed- 
dotted line is the conductance in first order in the tunnel 
coupling. As it is well known, the appearing peaks cor- 
respond to the addition energies of the dot being at reso- 
nance with the chemical potential of the leads. The blue, 
full line is the full first-order plus second-order result. For 
comparison, we also show (red, dashed line) the sum of 
sequential- and cotunneling under the neglect of renor- 
malization corrections. In Fig. [3{b) only the separate 
second-order corrections are displayed. As expected, we 
observe from these figures that only cotunneling effects 
lead to finite contributions in the Coulomb-blockaded re- 
gions while renormalization corrections contribute to the 
conductance only close to the resonances. 

Suppose, one wants to determine the renormalizations 
from the analytical expression of the linear conductance 
calculated from the diagrammatic real-time approach di- 
rectly. Once we have written the result in the form 
of Eq. (|3ip . we can read off Se, 5T ate and STa^+u- 
This is, however, not the only possibility. Supposing an 
overall renormalization of the tunnel coupling strength, 
without taking account for the different renormaliza- 
tions 5T a e and <5r Qje+ ;y of the two resonances, quantum- 
fluctuation effects on the tunnel coupling were described 
by one tunnel-coupling renormalization ST a only, i.e., 

(2) - 9% + + £ a ^5T a . This description 
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de uc T dr 

is indeed possible and leads to 



l-f(e)-f(e + U) 



5T L 



(33) 



l-f(e) + f(e + U) 

with ST a ^ from Eq. (f28|) . The above identity is eas- 

^ (i) 

ily understood by combining the condition ^ r ST a = 

— ST at€ + jjp- 2 — -8T a ^ + u with the relation 8T a ^ e = 
— ST a ^ e+ u that we found via the canonical transfor- 
mation. The expression for the single tunnel-coupling 
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FIG. 3. (color online) (a) Dimensionless conductance as a 
function of the level position in units of F in first order (black, 
dashed-dotted line) , including second order corrections due to 
real cotunneling (red, dashed line), and including all second 
order corrections (blue, full line), (b) Separate second-order 
corrections: real cotunelling (pink, dashed-dotted line), level 
renormalization (green, dashed line) and r-renormalization 
(violet, full line). We choose IY = Ir, k B T = 2T and U = 

3or. 



renormalization Eq. (|3"3")l is plotted in Fig. 2] for the 
zero-bias case; this total result is electron-hole symmet- 
ric. Also in this case, the ratio 8T a /T a is the same 
for both leads and we drop the lead index in the fig- 
ure. Far away from the particle-hole symmetry point, 
e + U/2 > k B T or e + U/2 < -k B T, only one of the 
transitions between the empty and the single occupation 
or between the single and the double occupation plays a 
role. As a consequence, ST a is given by the corresponding 
tunnel-coupling renormalization, <5r a . £ or 5T a . e +u, only. 
Around the particle-hole symmetry point, <5r Q interpo- 
lates between ST a ^ e and ST ate+ ij by averaging with the 
proper weights. 

In Refs.[28| and[2j| we calculated the relaxation rates of 
a quantum dot brought out of equilibrium, performing a 
perturbation expansion within a real-time diagrammatic 
approach. Similar as for the linear conductance discussed 
above, the second-order corrections could be fully under- 
stood in terms of cotunneling processes and renormaliza- 
tion of energy and tunnel coupling. Ignoring the possi- 
bility of different renormalizations 5r Qj£ and 8T ate +u, we 
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FIG. 4. Combined renormalization of the line width ST as 
function of the level position e in units of F, for zero bias 
voltage. We choose k B T = 2T and U = 30I\ 



extracted from the relaxation rates an expression for ST a 
that coincides with the one obtained from the differential 
conductance given in Eq. (|33l) . 



V. CONCLUSION 

We present an approach to classify the effects of 
quantum fluctuations in quantum-dot systems within a 
perturbative expansion in the tunnel-coupling strength. 
This approach is based on a canonical transformation 
that removes off-shell parts of the Hamiltonian and, si- 
multaneously, generates new transitions as well as renor- 
malizes system parameters such as energy and tunnel 
coupling. We illustrate this idea for the example of a 
single-level Anderson impurity model. Most notably, we 
find that the tunnel coupling strength for the two reso- 
nances connecting empty and single occupation respec- 
tively single and double occupation of the quantum dot 
renormalizc with opposite sign. The discussed effects are 
identified in a full second order expression for the linear 
conductance through the interacting quantum dot. 
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C. Energy renormalization in first-order transport 

For the simple single-level Anderson model considered 
in this paper, the energy renormalization gives rise to 
corrections in second-order transport but does not show 
up in the conductance in first order in V. The reason 
is that in the quantity of interest, Eq. only the 

equilibrium probabilities to zcroth order in T enter: the 
corrections linear in V happen to drop out. This is, how- 
ever, not always the case. For slightly more complex sys- 
tems, tunncling-induccd energy renormalizations already 
affect lowest-order transport. In a single-level quantum 
dot attached to noncollinearly magnetized ferromagnetic 
leads, e.g., there is a spin-dependent energy renormaliza- 
tion that has been described as an effective, tunncling- 
induced exchange field that influences the dynamics of 
the quantum-dot spin in a way that is detectable already 
in the linear conductance to lowest order in the tunnel 
coupling^ This happens also for a metallic island cou- 
pled to noncollinearly magnetized ferromagnetic leads i^I 
A similar effect is found for the case of double quantum 
dots ; 14 ' 24 ' 31 in which the two levels in the two dots de- 
fine an isospin that experiences a similar pscudo exchange 
field. The effect of effective exchange fields due to level 
renormalization was also observed in carbon nanotubes 
with orbital-dependent tunnel couplings 3 ^ and molecular 
single-electron transistors^ 3 - Furthermore, the first-order 
transport characteristics through quantum dots attached 
to superconducting leads^ shows features of different 
tunncling-induccd energy renormalizations of the empty 
and the doubly occupied dot. 



Appendix A: Cotunneling: Analytic results 

In this section we give the expressions for the Fermi 
golden rule rates in terms of digamma functions and their 
derivatives. We find for the rates describing a transition 
leaving the state of the quantum dot unchanged 
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The spin-flip term is evaluated as 



r v 

1 a L a 



1 



[* a (e,C0-*a'(e,E0] 



with $ Q (e, U) defined as in Eq. (|28[) . Finally we give 
explicit expressions for the cotunneling rates where the 
state of the dot changes by a charge 2e, 



/ F F / 



[$ a (e, U) + $ a , (e,U)] 
[$ a (e,U) + ® a ,(e,U)} 



All cotunneling terms can hence be expressed via the 
same digamma functions and their derivatives that ap- 
pear in the renormalization contributions. In order to 
calculate the current or the conductance, these rates have 
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to be weighted with the respective factor taking account for the number and direction of transferred particles and 

with the probability that the initial state is occupied. 
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